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Abstract 

In this paper we calculate at high-precision the expansions in e = (4 — D)/2 of the master 
integrals of 4-loop vacuum bubble diagrams with equal masses, using a method based on the 
solution of systems of difference equations. We also show that the analytical expression of a 
related on-shell 3-loop self-mass master integral contains new transcendental constants made 
up of complete elliptic integrals of first and second kind. 
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One popular approach to calculating Feynman diagrams relies on the use of integration-by-parts 
identities[l, 2]. By using algebraic identities the contribution of a diagram is reduced to a combination 
of irreducible integrals, the so-called master integrals. Only for these an explicit calculation is 
needed. For diagrams containing different masses and momenta, it is possible to express a multi- 
scale master integral in terms of single-scale integrals by means of well-known technique of asymptotic 
expansions[3, 4, 5, 6, 7]. This technique has had numerous applications in the last years; we refer 
the reader to the reviews [8, 9, 10] and the references therein. 

The single-scale master integrals contain at most two possible values of masses, and to, and 
external lines on-shell; being pure numbers, they can be calculated analytically or numerically at 
high precision once and for all. Some results are already available in the literature: two-loop single- 
scale self-mass master integrals and three-loop vacuum master integrals were evaluated analytically 
respectively in Refs.[12, 13, 14, 15] and Ref.[16], for all possible combinations of masses and to; 
among the numerous applications, for example, three-loop vacuum integrals were used for calculating 
the relation [11] between the MS and the on-shell quark mass at order a^. 

The three-loop single-scale self-mass master integrals appearing in the 3-loop QED contribution 
to the electron g-2 were calculated in analytical form in Ref. [17, 18]; other three-loop self-mass 
master integrals appearing in QCD diagrams were calculated analytically in Ref. [19]. An analytical 
calculation of self-mass master integrals which appear in the 4-loop QED contribution to the electron 
g-2 seemed not to be feasible with the techniques used so far, so that the author has recently developed 
a new method of calculation of master integrals based on the high-precision numerical solution of 
systems of difference equations [20, 21]. As first applications of this method, the three-loop all- 
massive self-mass master integrals and deeper expansions of 3-loop QED master integrals have been 
calculated at high precision [20, 21, 22]. 

In this paper we apply this method for the first time at 4-loop level, by calculating all the 4-loop 
vacuum-bubble master integrals with all massive lines. We think the results obtained are of some 
importance in view of the applications to 4-loop g-2, since the values of the these vacuum integrals 
can be used as initial conditions for the integration of differential equations in masses [23] and external 
momenta [24] in order to work out the 4-loop g-2 self- mass master integrals. 

We consider here only the master integrals which do not factorize into a product of integrals with 
fewer loops, already known. There are 10 different topologies for these integrals, shown in Fig.l, 
of which that corresponding to V13 is not planar. Three diagrams have two master integrals. The 
Laurent expansion in e = (4 — D)/2 of the master integrals are: 

V 1 = -2.5e" 4 - 11.6666666666666666666666666667e" 3 

- 31. 7013888888888888888888888889c" 2 - 67.5289351851851851851851851852e" 1 

- 140.2205432875405077693377105 - 573.5347004606566057988634367e 

- 2756.21982203281444579061754077e 2 - 18239.9256745938582375474975495e 3 

- 86167.478580009225683821281651e 4 - 468163. 766003582350760458083774e 5 

- 1976611. 1 12244975471 16753862679e 6 - 9573236.83728658266329080404209e 7 + 0(e 8 ) , (1) 
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V 2 = -1.6875e~ 4 - 7.8125c' 3 - 21.2096354166666666666666666667e" 2 

- 44.7695529513888888888888888889^ - 97.0765348083978788434432266577 

- 290.923474358239921908941634044c - 1719.80956172642828117680003592e 2 

- 8934.7307852874907107364237238e 3 - 51529.168884049820479322933243k 4 

- 236942.747528861416676127442283e 5 - 1 169551. 39766202377654297049171e 6 

- 4923196.24852687778836136788114c 7 + 0(e 8 ) (2) 

V 3 = 1.625e~ 4 + 10.2291666666666666666666666667e" 3 

+ 33.4657920462820508294798181834e~ 2 + 83.6384721220072849057411859385c- 1 
+ 142.482921255630165397059484622 + 73.8226346988098680135859408617e 

- 2215.06968731839973043724485006e 2 - 10483.810544518473050561125786c 3 

- 75868.174172897651897053047550k 4 - 266235.221240550768248796525565e 5 

- 1562163.76406718387868426308061e 6 - 5089761.51470684069337522895124c 7 + 0(e 8 ) , (3) 

V 4 = 1.5c" 4 + 9.5e~ 3 + 33.5e~ 2 + 59.8938292905212171438007855155c- 1 

- 6.77093494521947617804741496782 - 927.184486122109949342449134087c 

- 7607.58935727331640672412576731e 2 - 36043.9315389277686247446140344e 3 

- 199692.43129042192811514839511e 4 - 784516.116603849522465732718633c 5 

- 3938768.43031938169221908816973e 6 - 14490248.5601255703505458679593e 7 + 0(e 8 ) , (4) 

V 5 = -0.25e~ 4 - 2.33333333333333333333333333333e~ 3 

- 11.7695293784587392828307280615c" 2 - 47.315515044289411703624686853k- 1 

- 139.695412330272368658893294367 - 952.24968928016088 14725513796k 

- 1559.77276310374672922449535587e 2 - 17433.275615591257971452720325k 3 

- 18039.431294353693419505748424e 4 - 300299.661276324759969342746297e 5 

- 223362.891943298104957627580543e 6 - 5002799.01434332803606709511202c 7 + 0(e 8 ) , (5) 

V Q = 2.375e~ 4 + 15.9791666666666666666666666667e" 3 

+ 55.370266375339080107172526045e- 2 + 180.944159673951583547109549717c- 1 

+ 496.624640289958679240601686338 + 2103.52942297240529902308231413e 

+ 4385.91019392343391919235019438e 2 + 29880.7052368668975066622604109e 3 

+ 45750.0902911026047573732270476c 4 + 472601.136971716877928457965404c 5 

+ 573888.829426582538928155340449c 6 + 7686577.54749127119979436449377c 7 + 0(e 8 ) , (6) 

V 7 = -0.5c" 4 - 4.41666666666666666666666666667e- 3 

- 18.4682782839316233281953131334e~ 2 - 68.0334660202028428854058007578c" 1 

- 295.653297261799619866589394558 - 1111.65699468499496145662198062c 

- 4624.23997399795582698262827805c 2 - 18166.6470161776485664443287094c 3 

- 72747.4433773622788330319854055c 4 - 294544.5996338766980685312612k 5 

- 1152092.06523577988360701658787c 6 - 4748621.41279753209699735445982c 7 + 0(e 8 ) , (7) 
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V 8 = -0.75e" 4 - 6.25e~ 3 - 23.7024174258974349922929697001e" 2 

- 83.9541325752545325454615616904c" 1 - 438.152668276175094650239945304 

- 1151. 7188201038655941347197846c - 7624.47573539137683288770528724e 2 

- 16496.7827202671608148012112968c 3 - 127976.07440527009995741832362c 4 

- 245588.151837894378828676455499c 5 - 2103545.82723681108089748244607c 6 

- 3759879.03947407762892613884153c 7 + 0(e 8 ) , (8) 

Vg = 1. 80308535473939142809960724227c" 2 - 9.0439355503127513874520962153c" 1 

+ 42.0771149080147860090872289595 - 185.506435563765087317986949292c 

+ 794.690451268783109449998200836c 2 - 3340.85301017420853074405241137c 3 

+ 13860.9560593967975329182060014c 4 - 56963.5538585614643909090940607c 5 

+ 232474.430355572470669990540268c 6 - 943890.666478469607638755776088c 7 + 0(e 8 ) , (9) 

V w = 5.18463877571684963165682743229c" 1 - 38.9591108409575465929423742741 
+ 216.540829715907801334069893492c - 1056.20279561819887173122551098c 2 
+ 4804.63905121999951090192731115c 3 - 20965.3713155251719483741628451c 4 
+ 89093.2196878747322139561202477c 5 - 372004.961361167636317053811767c 6 

+ 1534700.00918085457184730424284c 7 + 0(e 8 ) , (10) 

V u = -1.36872356982626095206640482818c" 2 - 10.1337693663646787142726963524c" 1 
+ 25.6153496410310292617132790055 - 334.538889614391947948784653628c 
+ 984.338776062634546439769565907c 2 - 7183.20521761939021795215293201c 3 
+ 21442.9895505509641128293112066c 4 - 132391.959284919922968209921626c 5 

+ 396184.588583865634141091817143c 6 - 2278903.80544414662394641296011c 7 + 0(e 8 ) , 

(11) 

V 12 = 1.34894802170970895986445430292 - 8.34992524160927013803891299755e 

+ 40.4425656120211378353957075462e 2 - 177.294826919557277843981859823c 3 

+ 742.428926688827707784575315859c 4 - 3039.49399009410599098768805594c 5 

+ 12302.4392626409424372470869128c 6 - 49506.0256176923061393145333901c 7 + 0(e 8 ) , 

(12) 

V 13 = 0.997672576874263051049093586736- 5.83175720330362704111639113258e 
+ 27.5586906484187171474221359635e 2 - 119.248620281553405809848697481c 3 
+ 495.915629767162238727475039832c 4 - 2022.83520871375744081923368263c 5 
+ 8171.6838591733877817993555444c 6 - 32850.412079699027325130879249c 7 + 0(e 8 ) . (13) 

The master integrals Vj are defined as 

y. = (n D / 2 T(l + e))~ 4 / d D h d D k 2 d D k 3 d D k 4 , 
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where P 2 = (hi ■ k 2 ) 2 , P& = (fci + k 2 ) 2 , Pn = k\ ■ k 2 and all other Pj = 1; Qj is the product 
of denominators of the corresponding jth diagram of Fig.l, each with unit mass. The expressions 
of the finite integrals V\ 2 and V13 have been checked by comparing them with values obtained by 
performing Monte-Carlo integrations over Feynman parameters. 

The rapid growth of the coefficients is due to the fact that the integrals, seen as functions of D, 
have poles near D — 4. For example, in the case of V13 the nearest poles are D = 9/2, 5, 16/3, and 
11/2; by factorizing out these poles, one obtains a series with decreasing coefficients: 

(^.9976726 + 2.981017e + 2.981995e 2 
+ 1.315019e 3 + .432643e 4 + .0999339c 5 + .0150665c 6 + .0086152c 7 + 0(e 8 )^ . (14) 

Now we sketch here the method used for obtaining Eqs.(l)-(13); for a complete description the 
reader is referred to Refs.[20, 21]. For each master integral Vj we have chosen a denominator and 
raised it to n, and we have built a system of difference equations between the functions Vj(n). The 
difference equation for a given master integral contains in the r.h.s. only master integrals with less 
denominators, which are simpler. The equations of the system are solved one at once, beginning 
with that corresponding to the simplest master integral, Vi, and ending with that corresponding 
to the most complex, V13. The necessary boundary conditions at n — > 00 are expressed by 3-loop 
integrals which are calculated by solving other systems of difference equations. The homogeneous and 
nonhomogeneous solutions of the difference equations are expanded in factorial series or alternatively 
transformed into integral representations, and evaluated numerically for suitable integer values of 
the exponents. By taking these results as input values and by using the equations of the system as 
recurrence relations, one recovers the values of the master integrals V}(1), Eqs.(l)-(13). 

For example, the function Vi(n) satisfies the fourth-order difference equation 

p 1 V 1 (n+l)+p 2 V 1 (n)+p 3 V 1 (n-l)+p A V l (n-2)+p 5 V l (n-3) = Q(D - 2) 3 J 3 (l) J(n - 2) , (15) 

Pl = - 90n{n - l)(n - 2) , p 2 = 9(n - l)(n - 2)(16n - 9D - 6) , 
p 3 ={n - 2) (-20n 2 + (122 - 2lD)n + 27 D 2 - QQD - 42) , 

p 4 = - 32n 3 + (54 + 93D)n 2 + (-78L> 2 - 153L> + U)n + 18/J 3 + 87L> 2 + 21D - 24 , 
Pb = -(n-D- l)(n -2D + l)(2n - 3D) . 

It contains in the r.h.s. the integral obtained from V\ by contracting a line, which factorizes into 
a product of 4 one-loop tadpoles. The solution of Eq.(15) compatible with the large-n boundary 
condition V\{n) oc n~ D l 2 can be written as V\{n) = C I HO (n) + I NH (n) , where I HO and I NH are 
respectively the solution of the homogeneous equation with above large-n behaviour and a particular 
solution of the nonhomogeneous equation (15). The constant C turns out to be the value of the 
3-loop vacuum bubble obtained from V\ by deleting a line. 

The amount of calculations needed to work out and solve all the systems of difference equations 
is rather high, so that the calculations have been performed by means of an automatic tool, the 
program SYS described in Ref.[20]. 

Eqs.(l)-(13) have a precision sufficient to substitute analytical expressions for most of practical 
uses. Anyway, it is interesting to work out the analytical expressions of some of these results. In the 
case of V4, a straightforward application of the methods used in [17] for calculating the 3-loop g-2 



V 13 = 



(l+4e)(l + 2e) 



'4 



>4 
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master integrals gives 



V, = \e-± + ^ + ^ + (-3C(3) + ^) - ^ + 97C(3) - ^ + e(l89 C (5) 

_ _ 32vr 2 i n 2 2 + 768Li 4 (I) + 32 In 4 2 + 1261C(3) - ^) + 0(e 2 ) , (16) 

OU \ Z J Z J 

which agrees perfectly with Eq.(4). As usual, the expression turns out to contain integer values 
of Riemann ^-function ((p) = YlnLi 1 / nP > anc ^ vames of the polylogarithms Li p (x) = Y^=i xn l nP ■ 
The same method may be also used for calculating Vg, but not for calculating the remaining master 
integrals. The reason is the appearance of double elliptic integrals which cannot be expressed in 
terms of polylogarithms. 

The appearance of elliptic integrals in the calculation of multiloop diagrams with many massive 
lines is a well-known phenomenon. In fact, the principal difficulty of the analytical calculation of 3- 
loop g-2 has been the individuation of the approaches which allowed us to avoid elliptic integrations. 
Unfortunately in the case of the 4-loop vacuum master integrals the problem cannot be circumvented 
and we are forced to perform the elliptic integrations. 

Let us consider the simplest vacuum diagram, V\. By cutting one line of this bubble one obtains 
the 3-loop self-mass diagram shown in Fig.2. If we take all masses equal to one, and the external 
line on the mass-shell, this diagram has two master integrals, whose S\ is 

S x = (vr^l + e ))~ 3 J d D h d D k 2 d D k 3 I , 

and Q is the product of the four denominators of the diagram. 

The analytical expressions of V\ and Si are strictly related. In fact, introducing hyperspherical 
variables [25, 26] and performing the angular integrations one finds that the finite parts of V\ and Si 
can be expressed as a sum of many double elliptic integrals with the same three square roots in the 
denominator. The simplest of these integrals is that with unity as numerator 

oo oo 

A— f ——^ [ — ^ = 2.641379476074689 . . . , (17) 

J R(l, -1,-1) J R(m,l,-l)R(m,-l,-l) v ; 

o o 



R(x,y,z) = \J x 2 + y 2 + z 2 — 2xy — 2xz — 2yz . (18) 

We have been not able to calculate this integral in analytical form. Therefore, we have resorted to 
'experimental mathematics': we have evaluated it at very high precision and we have tried to fit the 
numerical value with various kinds of analytical expressions. Luckily 1 , we have found that 

A = K(w-)K(w+) , W ± = ^TT> z± = -{2-V3)\4±Vl5) 2 , (19) 

z± l 

1 We have been inspired by the analytical expressions of a double elliptic intcgral[27] related to the Appell F 2 
function and by the analytical result of the 2-loop 3-lines "sunrise" self-mass diagram in two dimensions 

S(p 2 = -1, ml = 1, m\ = 1 - 2- 1 / 2 , m\ = 1 + D = 2) = itK(1/2)/2 = V^r 2 (l/4)/8 . 
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where K is the first of the two standard elliptic integrals 



i 



1 



K(m) = 



dt 




dt VI - mt 2 



(20) 



VT^Vl - mt 2 ' 



o 







We have verified Eq.(19) with a precision of more than 30000 digits 2 . Once the analytical expression 
of the basic integral A has been identified beyond any reasonable doubt, we have considered the 
full integral representations of the finite parts of V\ and Si; they contain integrals which differ from 
Eq.(17) by additional polynomials and logarithmic functions in the numerators. A high-precision 
numerical expansion in e of S\ was calculated in Ref. [28] by means of the program SYS. We have 
been able to fit the numerical result with the analytical expression 



by using the integer-relation search algorithm PSLQ[29]. We have verified Eq.(21) with 1200 digits of 
precision. Eq.(21) is important because it shows for the first time that there exist single-scale integrals 
which contain non-polylogarithmic transcendental constants. The same constants also appear in the 
analytical expression of V±, together with higher-transcendentality constants so far unidentified. 

This result also sheds light on the possible analytical structure of of the 4-loop QED contribution 
to the electron g-2. In fact one of the 4-loop g-2 master integrals is a 4-loop 5-lines self-mass integral 
analogous to Si, expressible in terms of irreducible triple elliptic integrals analogous to Eq.(17). 
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Figure Captions 



Figure 1: Four-loop vacuum bubble diagrams. 
Figure 2: Three-loop self-mass diagram. 
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